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PROBABILITY  DISTRIBUTION  OF  ARRAY  RESPONSE 
FOR  RANDOMLY  PERTURBED  ELEMENT  GAINS 

INTRODUCTION 


Although  an  array  beamformer  with  fixed  element  positions  may  be  de¬ 
signed  for  good  side  lobes  or  skirt  selectivity  behavior  by  choice  of  the 
element  weights,  the  actual  array  will  undergo  degradations  due  to,  e.g., 
element  position  movement,  element  gain  quantization,  random  element  gains, 
etc.  Here  we  will  investigate  the  effect  of  randomly  perturbed  element  gains, 
each  with  omnidirectional  response,  on  the  power  response  of  the  array  beam- 
former  for  a  single-frequency  plane-wave  arrival.  (We  do  not  consider  the 
effects  of  random  phase  shifts.)  Since  the  power  response  to  a  single-frequency 
plane  wave  is  itself  then  a  random  variable  at  each  angle  of  look,  we  will 
evaluate  its  mean  and  variance  as  a  function  of  the  element  locations,  the 
look  angle,  and  the  plane-wave  arrival  frequency,  propagation  speed,  and 
direction.  More  generally,  we  will  evaluate  the  characteristic  function  of 
the  random  power  response.*  From  these  results,  we  will  be  able  to  make  quant 
tative  statements  as  to  the  required  tolerance  on  the  individual  element  gains 
for  a  specified  side  lobe  level,  as  a  function  of  the  number  of  elements  and 
their  weight  structure. 


GENERAL  ARRAY  CONSIDERATIONS 

The  geometry  to  be  utilized  here  is  illustrated  in  figure  1. 


Z 


y 


Figure  1.  Geometry 


-it  < 


♦Although  not  pursued  here,  a  single  numerical  Fourier  transformation 
would  then  yield  the  probability  distribution  of  the  power  response  at  the 
particular  look  direction,  arrival  angle,  etc.,  of  interest. 
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The  radial  vector  to  the  location  of  the  k-th  element  of  an  arbitrary  array 
(numbered  in  any  convenient  fashion) is 

■vH  *  *  ^7  +  **  k 

=  rK  si**4  (cose^T  +  5)n  ~$)  +■  ^  CU  k  , 

(1) 


where  “C,  J  J(  are  unit  vectors  in  the  directors  of  the  x,  y,  z  axes,  respec¬ 
tively.  The  arrival  direction  vector  of  a  plane  wave  from  direction  (£,^) 
is 


U  *  3ivi 


(cos  0  V  v  0  j)  +  ««  cj>  k 


(2) 


If  we  adopt  the  convention  that  the  time  delay  is  zero  for  the  plane 
wave  to  reach  the  origin  in  figure  1,  then  the  time  delay  to  the  k-th  element 
is  (for  speed  of  propagation  c) 


-  iu v*  -  -  £-[ykcose  j'<« ^ + y* si* e  ?>* <*s $). 


(3) 


Also,  if  the  look  direction  of  the  array  beamformer  is  v*k  ,  then  the  delay 
employed  in  the  k-th  receiver  element  (to  accomplish  maximum  response  when 
steered  in  the  arrival  direction)  is 

3m4  +  *«cos  4] .  (4) 

If  the  plane  wave  arrival  is  characterized  by  the  waveform  s(t)  at  the  origin 
in  figure  1,  then  the  waveform  at  the  output  of  the  delay  element  in  the  k-th 
receiver  channel  is 

5  [i  4-  4)]  .  (5) 

For  a  beamformer  that  weights  this  k-th  waveform  by  w*  and  sums  over 
the  available  element  outputs,  the  voltage  transfer  function  at  frequency  f 
of  the  beamformer  (as  applied  to  the  plane-wave  waveform  s(+))  is 

A  [t ,  »M,  i')  *  «p  [- '  fe  ,*)  -  ■ T« ^ €] .  161 


2 


■  .  ..  -  . 
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This  is  a  general  expression  for  the  voltage  transfer  function  of  the  beam- 
former  of  an  arbitrary  volume  array.  We  will  adopt  an  abbreviated  version 
of  (6)  for  later  notational  convenience: 


A-  -  » 

where 

VK  =  \  eXj>  [j  I  2-rr-f  ~ 

is  a  deterministic  complex  function  of  angles  ,  frequency  f,  and 

element  positions  • 

LINEAR  ARRAY 

The  transfer  function  in  (6)  applies  for  an  arbitrary  array.  In  the 
special  case  of  a  linear  array,  we  set  (see  figure  1) 

\  -  %  -  °> 

thereby  obtaining 

c10) 


(7) 


(8) 


and 

^  *  WK  »  2irf -§<-  (cos  ^  -  COS  <j>)]  .  (11) 

Broadside  arrivals  correspond  to  $>  =  tt / 2 . 

For  the  further  specialization  of  an  equi-spaced  line  array  of  spacing  d, 
we  let 


2k  =  K  J, 


obtaining 


VH  e*j.(-ik«), 


where  dimensionless  variable 


U  *  2-rrf  ^  (cos  -  cos  (£)  . 


(12) 


(13) 


(14) 
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Then  voltage  transfer  function  (7),  for  this  special  case,  becomes 

A  =  ^Wke*?(-iku).  (15) 

Since  exp(-iku)  has  period  2tt,  there  is  no  need  to  investigate  A  outside 
the  range  (-tt,  tt)  of  u .  This  range  will  cover  all  possible  cases  of  element 
spacing,  look  direction,  plane-wave  arrival  angle,  frequency,  and  speed. 

PLANAR  ARRAY 

In  the  special  case  of  a  planar  array,  we  place  the  array  in  the  x-y 
plane  by  setting  z^  =  0;  see  figure  1.  Also,  we  limit  consideration  to  a 
planar  array  with  a  grid  structure;  see  figure  2.  That  is. 


Figure  2.  Grid  Structure  of  Planar  Array 

at  each  and  every  x-coordinate  xm,  the  y-coordinates  of  the  element  positions 
are  {y  };  the  spacing  need  not  be  uniform  in  either  coordinate  direction. 
Using  a  more  appropriate  labeling  system,  we  express  (7)  as 


where  (from  (3)  and  (8)) 


Tmr>(£,  <|>)  =  COS  e  +  Vj„3in  ©]  (17) 

id 

=  2irf ^063^-  Sl»4  toS0)  4-  -^St*  ^  St*  -  S'"  ^  SIM $)] J  (18) 
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For  the  further  specialization  where  the  weights  are  obtained  via  a 
multiplicative  pattern, 


*Nl'\ 


(19) 


and  (16) -(18)  yield  the  voltage  transfer  function  as 

l\  ,  -T-(5m<ljCCJ  6^- cose)]  • 

£ i 2rf 4  (SH  5'”  §  -  S’"  +  s'"  6)]  *  i'1-  ?v-’  •  (J0) 

Thus  the  transfer  function  in  this  case  can  be  represented  as  the  product  of 
two  component  responses.  If  weights  {w  M }  are  chosen  such  that  the  first 
component  in  (20)  has  good  side  lobe  behavior*  and  weights  {w^  }  are  similarly 

chosen  for  the  second  component  in  (20),  then  response  A  will  have  good 
side  lobe  behavior  for  all  look  and  arrival  angles. 


Finally,  when  the  elements  are  equally  spaced  by  dx  in  the  x-direction 
and  by  dy  in  the  y-direction,  (20)  yields  for  the  voltage  transfer  function. 


A 


(21) 


where  dimensionless  variables 

U.  »  2-rrf  ^ [s'm  ^  Cos  6^  -  SIM  <j>  cos 

v  ^  2-rr-f  (s\*  sin  ^  51*1  $). 


(22) 


Since  (21)  has  period  2tt  in  u  and  v,  there  is  no  need  to  investigate  A  outside 
the  range  (-7T,  tt)  in  u  and  in  v. 


STATISTICS  OF  POWER  RESPONSE  OF  ARRAY 
FOR  GAUSSIAN  PERTURBATIONS 


The  general  voltage  transfer  function  was  given  in  (7)  and  (8) ,  where 
{w^}  are  the  design  weights  of  the  array  beamformer  and  are  assumed  real. 
However,  because  of  imperfections  in  manufacture  or  control,  the  gains  of 
individual  receiving  elements  are  not  equal,  but  are  randomly  perturbed  about 
their  design  values.  We  incorporate  this  feature  in  our  voltage  transfer  func¬ 
tion  (7)  by  saying  that  the  modified  (perturbed)  voltage  transfer  function  is 
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A  -£(■  +  >!)*,  (25) 

where  is  a  (dimensionless)  real  random  variable  (with  zero  mean  without 
loss  of  generality)  that  measures  the  relative  perturbation  of  the  k-th 
element  voltage  gain.  We  will  assume  that  {r^}  are  joint  Gaussian  with 
correlation 

^{^^3  =  »  (24) 


that  is,  we  allow  different  elements  to  have  correlated  gain  perturbations. 
(Later,  the  Gaussian  assumption  will  be  eliminated  in  adless  complete  statis¬ 
tical  analysis.  Also,  even  if  (r^l  are  not  Gaussian,  A  will  tend  toward  a 
complex  Gaussian  random  variable  for  large  numbers  of  elements,  by  the  Central 
Limit  Theorem.  Thus  the  following  results  have  applicability  for  a  wider  range 
of  cases  than  originally  presumed.) 

The  power  transfer  function  of  the  perturbed  array  is 


p  S|sr-|2(i^)' 


(25) 


We  are  interested  in  the  probability  distribution  of  this  real  random  variable. 
In  appendix  A,  the  distribution  is  shown  to  be  given  by  a  double  integral  of 
a  correlated  second-order  Gaussian  density  (with  unequal  variances)  over  an  off¬ 
set  circle.  A  limited  table  of  this  quantity  is  available  in  reference  1; 
however  it  is  not  extensive  enough  for  our  purposes  and  it  yields  little  in¬ 
sight  into  the  array  behavior.  Therefore,  instead,  we  will  derive  the 
characteristic  function  of  random  variable  p,  and  from  this  quantity 
extract  the  cumulants  of  P.  The  characteristic  function  of  p  is  derived 
in  appendix  B  and  is  given  by 

f(?)=  =  [l-i? 2D +(>?)’ (j*- |s)’j]  h ■ 


ex 


? 


I  -  >  f  2p  +  (» fT  (d’-  |Sl3)  J 


(26) 


where  A  is  given  by  (7)  and 

V  = 


%  v»  ■ 


S  =  2  2 V,  Vj  . 

*  J 


(27) 
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Thus  the  three  fundamental  functions  A,  D,  and  S  completely  characterize  the 
statistical  behavior  of  P.  A  single  numerical  Fourier  transform  of  (26) 
would  yield  the  probability  distribution  of  P  for  specified  values  of  A,  D, 
and  S,  if  desired;  see  reference  2. 

The  cumulants  of  P  are  derived  from  (26)  in  appendix  B;  the  first  three 
are  given  by 

E(f>)  -  IaT+  x>, 

=  I>3  +  |s(  +  2  dIA|’  4-  21fe{A15*]  ( 

'2®3+  fPlSf+fD'|Ar+  C|S|’lAr+  l2Dfc{A‘S*}.  (28) 


The  first  relationship  for  the  mean  E{P}  actually  holds  for  any  dis¬ 
tribution  of  {r^},  not  just  Gaussian;  this  is  most  easily  seen  by  direct  use 
of  (25),  (24),  (7),  and  (A-3) . 

Two  useful  special  cases  of  the  above  results  are  available  from  (7) 
and  (27): 


uncorrelated 

perturbations 


uncorrelated 
equal -variance 
perturbations 


(29) 


(30) 


k- 

d  =  <sw«r 

S  *  «ff»f  ^ 

where,  in  the  latter  case,  <^:  is  the  common  variance  of  the  random  gain  vari¬ 


able  introduced  in  (23) . 


STATISTICS  FOR  NON-GAUSSIAN  PERTURBATIONS 

In  this  section,  random  perturbations  {r^}  are  assumed  statistically 
independent  with  identical  statistics,  but  are  not  assumed  Gaussian.  We  let 


E{^}  £  tfl*  ;  Vn,  =  0 ,  =  <Tr  . 


(31) 
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The  mean  of  random  power  response  P  is  given,  by  use  of  (25),  (31),  and 
(7),  as 


£{P]  - 


^  J 


V  v. 
*  J 


*  |AP  +  *»  fKI  • 


(32) 


The  variance  of  P  is  derived  in  appendix  C  and  is  given  by 

V«r?  =  2  h,,  +  SefVSlvf}] 

+  4  w.  Tie  {  A  2.  Kf  V*] 

+  +  |Z-Vf]  +  kf. 

(33) 


2 

For  Gaussian  perturbations,  m3  =  0,  114  =  3m2,  and  (33)  reduces  to  (28), 
where  D  and  S  are  given  by  (30) .  For  a  uniform  probability  density  function 
for  {r,}  ,  we  have 

M,  =  0,  »|,=  Uir;,  (uniform) ,  (34A.) 

whereas,  for  a  zero-mean  exponential  density,  we  have 


(exponential) .  (34B) 


Generally,  when  m3  is  zero  (as,  for  example,  for  a  symmetric  density  for 
{  rfc  }),  the  only  change  in  the  mean  and  variance  of  P  from  the  Gaussian  re¬ 
sults  in  (28)  is  the  addition  of  the  1114  -  3m^  constant  term  in  (33);  the 
latter  two  examples,  above,  indicate  that  this  term  can  be  either  positive 
or  negative. 

One  specification  of  the  relative  randomness  r^  of  the  k-th  element 
voltage  gain  (introduced  in  (23))  is  its  variance  T*  .  However,  it  is 
sometimes  more  useful  to  express  this  variation  in  decibels.  This  problem 
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is  addressed  in  appendix  D,  and  results  are  plotted  in  figures  3  and  4.  For 
example,  from  figure  3,  for  Gaussian  perturbations,  if  ^  =  0.1  (<^x«0.0l  ), 
the  range  of  element  gain  perturbations  is  (-1.19,  1.05)  dB,  with  probability 
0.9  -  0.1  =  0.8;  that  is,  the  element  gain  is  less  than  -1.19  dB  from  design 
10  percent  of  the  time  and  greater  than  1.05  dB  from  design  10  percent  of 
the  time.  Figure  4  gives  results  for  a  uniform  distribution  of  element 
variations;  the  corresponding  range  for  <*  =  0.1  is  (-1.30,  1.13)  dB,  with 
probability  0.8. 


Figure  3.  Range  of  Element  Variations  for  Gaussian  Perturbations 
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Figure  4.  Range  of  Element  Variations  for  Uniformly 
Distributed  Perturbations 


APPLICATIONS  TO  ARRAYS 


When  the  perturbations  {r^}  in  (25)  are  not  zero,  the  random  pattern 
assumed  by  power  transfer  function  P  depends  on  the  particular  realization 
of  (rj, }  .  In  figure  5,  four  different  realizations,  for  a  20  element  equi 
spaced,  equi -weighted  line  array  with  ^  =  0.01,  are  plotted  for  u  in  the 
range  (0,tt  );  see  (12) -(14).  In  figure  6,  the  element  weighting  {wj<}  is 
changed  to  Hamming, 

%  =  &S+-  0.44  cos (2-rr^Tj-^  ,  IsfcsK, 

x 

and  the  variance  is  reduced  to  0.001.  Although  the  ideal  Hamming 
response  (r^  =  0)  would  continue  to  decay  with  u,  the  four  random  realiza¬ 
tions  in  figure  6  saturate  at  a  value  near  -40  dB  for  this  example. 


le  Array,  Equi-Weighted 
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In  figure  7,  four  random  power  responses  of  a  planar  array  of  16  x  16 
elements,  with  multiplicative  -30  dB  Doiph-Chebyshev  weighting  and  Of1* 
0.01,  are  plotted  for  v  =  0;  see  (19) -(22).  In  figure  8,  the  only  change 
is  to  set  v  =  u;  since,  according  to  (21),  the  ideal  pattern  drops  to  -60  dB 
along  this  si  ice  in  the  u,  v  plane ,  the  saturation  level  is  lower  in  this  figure 
than  in  figure  7. 

For  both  the  linear  and  planar  arrays,  we  wish  to  find  the  average  and 
the  standard  deviation  of  the  random  power  response  as  a  function  of  the 
number  of  elements,  weighting,  element  gain  variation,  etc. 


Figure  7.  Four  Random  Power  Responses  for  Planar  Array, 
Doiph-Chebyshev  Weighting;  v  =  0 
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Figure  8.  Four  Random  Power  Responses  for  Planar  Array, 
Dolph-Chebyshev  Weighting;  v  =  u 


LINEAR  ARRAY 


We  restrict  consideration  here  to  uncorrelated  and  equal -variance  gain 
perturbations,  an  equi-spaced  line  array,  and  symmetric  weights.  Then  (30) 
and  (13)-(14)  are  relevant  to  the  discussion.  We  find,  for  M  elements 
(M  even) , 


A 

3> 

S 


e,p(-i^u)A,(u), 

<■>’  Mo), 

CXf  5r(w), 


(36) 
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where  real  functions 

M"1)  •  2  2w.coi(ia|2^u), 

(J7) 

Then  (28)  yields,  for  the  statistics  of  the  random  power  pattern  P, 

ElP] «  A‘M  +  ^(o), 

2r;A;W[sv(oHS,W] 

+  «;4[sra(o)  +  5»], 

E{{f-  EfflT]  -  6  K  A’  (“)[^(°) +  ^ 

+  2<nfs.w[s;io)<-3s;(u0 . 

(38) 

Some  useful  properties  of  these  functions  are  listed  below: 

/\  -  exp^-imu)  -  Am(u)  =  Am(h+2-i»)  =  A„  (_v*), 

K  ^  *■  l 

M 

S  =  *  *~r  X<exj.(-«vn2u)5  -^(m)  =  5*(«+^, 

fy  ^  ^  Jv  Am  (v)  Am  (v-2u)  =  Si  (0, 

<T?  £  AM(V)  Am  (&"  v)  =  Si  M 


(39) 
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Thus  Aj^fu)  is  periodic  of  period  2tt  and  possesses  conjugate  symmetry;  SM(u) 
has  period  tt;  and  either  the  correlation  or  convolution  of  A^(u)  yields  SM(u) , 
where  2u  appears  as  the  shift  in  the  argument  of  the  correlation  or  convolu¬ 
tion.  These  properties  enable  us  to  limit  examination  of  (38)  to  the  range 
(0,7t)  in  u. 

In  figure  9,  an  M  =  20  element  line  array  with  equal  weighting  and  = 
0.01  is  considered.  Four  curves  are  drawn  in  the  figure;  starting  from  the 
bottom  up,  they  correspond  to  (1)  the  ideal  power  pattern,  (2)  the  mean  power 
pattern,  (3)  the  mean  plus  one  standard  deviation  (Std.  Dev.),  and  (4)  the 
mean  plus  two  standard  deviations.  Mathematically,  these  are 


(0 

(!' 

0) 

(+) 


£  ff-}  +  5fl  IVvlf} 


*ft>r  ^  ^  Q. 


(40) 


These  curves  are  labeled  m(  =  0) ,  m,  m  ♦  <r,  m  +  2  <r  ,  respectively.  Thus, 

for  u  near  tt,  although  the  mean  power  response  m  is  only  increased  by  about 
1  dB  (near  the  peak),  the  mean  plus  two  standard  deviation  response,  m  +  2<T  , 
is  increased  by  about  4  dB.  Probabilistic  interpretations  and  the  depth  of 
the  attainable  levels  in  the  deep  side  lobe  region  will  be  presented  later  in 
this  subsection. 
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In  figure  10,  results  for  a  20  element  array  with  Hanning  weighting  are 
presented;  that  is, 


I  <  M  . 


(41) 


Values  of  ^-=0.01,  0.001,  and  0.0001  are  given  in  figures  10a,  10b,  and 
10c,  respectively.  The  side  lobe  level  of  -30  dB  is  not  attainable  for 
=0.01;  whereas,  for  =  0.0001,  deep  side  lobe  levels  of  -45  dB  are 
attainable. 


Figure  10a.  Average  Power  Responses  for  Hanning  Weighting;  =  0.01 


0.001 


Figure  10.  Average  Power  Responses  for  Hanning  Weighting 
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Figure  10c.  Average  Power  Responses  for  Hanning  Weighting;  a '  =  0.0001 

Figure  10.  (Contfd)  Average  Power  Responses  for  Hanning  Weighting 


In  figure  11,  the  Hamming  weighting  of  (35]  is  employed.  It  is  seen 
in  figure  11a  that  =  0.001  gives  about  a  7  dB  increase  in  the  m  +  curve 
over  the  ideal  response  (at  the  peaks).  The  results  for  G?  =  0.0001  in 
figure  lib  yield  approximately  a  3  dB  increase. 


Figure  11a.  Average  Power  Responses  for  Hamming  Weighting;  =  0.001 
Figure  11.  Average  Power  Responses  for  Hamming  Weighting 
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Figure  lib. 


Average  Power  Responses  for  Hamming  Weighting; 


0.0001 


Figure  11.  (Cont'd)  Average  Power  Responses  for  Hamming  Weighting 


A  Dolph-Chebyshev  array  design  (reference  3)  is  considered  in  figure 
12.  This  particular  case  is  a  20  element  array  designed  for  -20  dB  side 
lobes.  Results  for  ^  =  0.01  and  0.001  are  presented  in  figures  12a  and 
12b,  respectively.  The  only  change  in  figure  13  is  that  the  design  is  for 
-30  dB  side  lobes;  it  is  seen  that  the  value  of  =  0.001  yields  a  2.5  dB 
increase  in  the  m  +  2«~  curve  over  the  m (<*>*=  0)  curve.  Finally,  in  figure 
14,  the  design  is  for  -40  dB  side  lobes,  and  <£*  =  0.0001  is  going  to  be 
necessary  to  take  advantage  of  this  design. 
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Figure  13a,  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  cr  =  0.01 


Figure  13b.  Average  Power  Responses 
-30  dB  Weighting;  a2 


for  Dolph-Chebyshev 
0.001 


Figure  13.  Average  Power  Responses  for  Dolph-Chebyshev  -30  dB  Weighting 


j 
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In  the  region  where  the  ideal  power  pattern  response  becomes  very  low  (co 
be  called  the  deep  side  lobe  region),  as,  for  example,  in  figure  10  for  larger 
u, or  in  any  of  the  figures  for  u  in  a  notch,  the  mean  response  curves  m, 
m  +  and  m  +  2<r  saturate  and  do  not  drop  continuously.  It  is  useful  to 

know  these  saturation  levels  and  what  they  depend  on,  so  that  design  criteria 
can  be  developed  for  future  cases.  From  (36)  and  (37),  we  have  (for  w  >0) 

=  1 

T>  =  <rr" 

K 

k  . 

**  —  \ 

(42) 

Therefore 

V 

Wl" 

(43) 

where 

>v  *  • 

s  M  \2 

(SO 

'  V»1*)  / 

±< 

M*  1 

(44) 

2 

Now  in  the  deep  side  lobe  region,  |  Ar«  D;  therefore  (28)  indicates  that  the 
depth  of  the  mean  power  response  is  approximately  D.  Since  the  curves  in 
figures  9-14  are  normalized  to  0  dB  at  their  peak  for  <Ty  =  0,  (43)  gives  the 
attainable  depth  of  t1  e  mean  power  response  relative  to  the  peak.  Table  1 
gives  the  values  of  ^e££  in  (44)  for  two  types  of  weighting.  By  Schwartz’s 
inequality,  Mef£  £  M;  the  upper  limit  is  attainable  only  for  equal  weighting. 


Table  1.  Values  of  M  rjC 

eff 


Weighting 

Meff 

Equal 

M 

Hanning 

M  *  1  fM  >  2) 
1.5 

i 

As  an  example,  for  equal  weighting,  M  =  20,  =  0.01,  then  (43)  and 

table  1  yield 


22 


JL  =  -2-^-  *  o.  boos'  -  -  33  <JB- 

Wl  20 


20 
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This  value  is  seen  to  be  attained  at  the  notches  in  figure  9.  A  second  ex¬ 
ample,  for  Hanning  weighting,  M  =  20,  ~  0*001,  yields 


2-  .  JL22L  ,  -  5  dB. 

IAC  2,/,y 


(46) 


This  value  can  be  seen  in  figure  10b  for  larger  u. 

2 

In  the  deep  side  lobe  region,  characterized  by  | A |  <<  D,  we  also  have 
| S | ^  <<  D2,  as  may  be  seen  from  (39)  and  (42) .  (This  is  true  except  for  u  near 
it.)  Then  (28)  yields 


E[p)  «  t>,  Var^f]  «  T>*. 


This  means  that 


Wt  * 

=  P  +  D  =  2'D> 
v»i+2<r  *  P  4-2T>*  3J>. 


(47) 


(48) 


That  is,  in  the  deep  side  lobe  region,  the  m  +<r  curve  is  3  dB  above  the  m 
curve,  and  the  m  +  2 <r  curve  is  4.8  dB  above  the  m  curve.  These  results  are 
independent  of  M,  <Tpa  ,  and  the  weighting.  The  behavior  predicted  by  (48) 
is  observed  in  figures  9-14. 

Equation  (43),  coupled  with  (48),  enables  a  simple  rule-of-thumb  for  the 
deep  side  lobe  region.  Furthermore,  a  simple  probabilistic  interpretation 
is  possible  in  that  region,  as  follows.  When  | A | ^  «  D  and  |s|  <<  D^,  (A-4) 

yields 


Then  (A-5)  and  (A-6)  yield 

TVoV.  (?<  0  21  1 '  e*f(-  b“)  >  L  *  °- 


(49) 


(50) 
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In  particular,  referring  to  (48) , 

<  M+ka)  =Prol.(iS<'J>(HV'j)  =  1-  exf(-l-k). 


(51) 


A  short  table  follows.  It  shows  that  the  probability  of  the  random  power 
response  being  less  than  its  mean  m  is  0.63,  less  than  m  +  <T  with  probability 
0.86,  and  less  than  m  +  2<r  95  percent  of  the  time. 


Table  2.  Probabilities  of  the  Random  Power 
Response  in  Deep  Side  Lobe  Region 


k 

0 

1 

2 

Prob 

0.63 

0.86 

0.95 

The  fact  that  the  probability  is  not  0.5  for  k  =  0  in  (51)  is  a  result 
of  the  asymmetry  of  the  probability  density  function  of  P  in  the  deep  side 
lobe  region*  in  fact,  (50)  yields  the  exponential  density  function  exp 
(-P/D)  for  P  >_  0.  (If  P  had  been  Gaussianly  distributed,  the  probabilities 
in  table  2  would  be  0.50,  0.84,  and  0.98,  respectively.) 

When  the  element  gain  perturbation  are  not  Gaussian,  (33)  indicates  that 
Var  {P}  is  changed.  However,  fo^  a  uniform  density,  the  results  in  (34)  give 
m3  =  0  and  1114  Proportional  to  <*V  •  For  the  small  values  of  of  interest 

in  the  array  examples  in  figures  9-14,  the  curves  for  the  uniform  density 
were  virtually  indistinguishable  and  have  not  been  presented. 


PLANAR  ARRAY 

For  a  planar  array,  the  random  gain  perturbations  change  (16)  to  read 


A  -  |^2(h-Ov..'  A+2Z 


r  v 

**»n  > 


(52) 


where  vmn  is  given  by  (18).  For  a  multiplicative  weight  pattern  (19),  (18) 
yields 


60  (a) 


(53) 


(xl  (vl 

where  vv  and  v  1  are  given  by  (20).  Thus 
m  n 
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For  an  equi-spaced  array  in  both  the  x  and  y  directions, 

w 


K  =  exp(-'iw,u), 

'Jf  *  €XJ»  (-1*1  v). 


(55) 


where  u  and  v  are  given  by  (22) . 


We  restrict  consideration  to  uncorrelated  and  equal-variance  gain  per¬ 
turbations.  Then  (30)  yields 

M  =  |?v»  1  , 

d  -  swrr  , 

Wi  n  "  * 

s-  -  o’  svf  xvf  .  <r;[^r'«r(-.2^)][^«P(-,n2v)l 


(56) 


For  the  particular  case  of  ^-9  =  0  ,  (22)  yields  v  =  0.  Equation  (56} 
then  reveals  results  that  are  identical  to  the  linear  arr^y,  except  that  | A | - 
is  scaled  by  r?*y  ,  and  D  and  S  are  scaled  by  ?*4’  That  is,  the 
line  array  results  are  directly  applicable  if  a  new  effective  element  toler¬ 
ance  of  /Nef f  is  used  insead  of  ,  where 


N 


zitr 


(57) 


Thus  the  extra  dimension  afforded  by  a  planar  array  causes  a  smoothing 
effect  on  the  random  power  response.  This  feature,  shown  here  only  for  a 
particular  slice  in  the  u,  v  plane,  will  be  shown  later  to  be  a  general 
result;  the  (generally)  increased  number  of  elements  in  a  planar  array  allows 
for  more  variation  in  the  element  gains. 

For  symmetric  weights  and  M  x  N  elements  (M,  N  even) ,  we  find 

A  =  exp(-  I  ^-vj  £  (u)  A?  (v), 

D  -  <r;  sr(0)5%\ 

5  =  «xf(-i(M+i)u-i^i)v)<r;  (Sg) 
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where  real  functions 

/Cm  ^  2  J\wcMp=^u) , 

|^/2  ^ 

5r%)  5  2  Cos((2*-M-l>), 


with  corresponding  definitions  for  real  functions  (v)  and  (v) .  Then 
(28)  yields  z  r  r 

EfP]  •  A ®(m)  A?  (v), 

Va4T)  -  2  rr2  Af  («) Af W[5T(o)Sr(o) + 5fw  srw] 

+<r;[s^(o)^‘(.)-sr'(u)^‘(o]. 


(60) 


fx) 

Some  useful  properties  of  these  functions  are  now  listed:  (u)  has 


(xT 


uuiiiv  vaoV'Xuj.  wi.  liivjv.  i.  uu^  Liuua  v  uun  i  xo  tvu  •  Aj*  y\uj. 

period  2tt  and  is  even  in  u;  A^yf (v)  has  period  2tt  and  is  even  in  v;  sj. 

has  period  2tt  and  is  even  in  u;  and  (v)  has  period  2tt  and  is  even  in  v 

Therefore  we  need  compute  (60)  only  fSr  0  <  u,  v  <  tt. 


(u) 


In  figure  15,  a  16  x  16  element  array  with  multiplicative  Dolph-Chebyshev 
weighting,  designed  for  -30  dB  side  lobes  in  both  x  and  y  directions,  is 
considered.  In  figure  15a,  =  0.1,  whereas  in  figure  15b,  <Tr  =  0.01; 

both  results  are  for  a  slice  along  the  v  =  0  axis.  It  is  seen  that  the 
element  gain  variability  should  be  at  least  as  small  as  that  in  figure  15b. 

In  figures  15c-15e,  the  slices  in  the  u,  v  plane  are  taken  along  the 
directions  v  =  77/8,  v  -  u/4,  and  v  =  u,  respectively .  The  -60  dB  level  in 
the  last  figure  for  larger  u  is  a  result  of  the  multiplicative  weight  patterns 
adopted.  The  average  responses  saturate  near  the  -40  dB  level  in  the  deep 
side  lobe  region. 


Results  for  an  array  of  equal  weights  are  given  in  figure  16,  whereas 
multiplicative  Hanning  weighting  is  considered  in  figure  17.  Both  results  are 
along  the  v  =  u  slice  in  the  u,  v  plane.  The  deep  side  lobes  predicted  by 
Hanning  weighting  are  filled  in,  as  figure  17  indicates.  The  pattern  for  a 
equi-weighted  array  is  not  affected  as  drastically. 


.*7* 


Figure  15a.  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array;  =  0.1,  v  =  0 


Figure  15.  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array 


Figure  15b.  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array;  =  0.01,  v  =  0 


Figure  15.  (Cont'd)  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array 


I 

Figure  15c.  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array;  o*  =  0.01,v  =  tt/8 


Figure  15.  (Contfd)  Average  Power  Responses  for  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array 
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Figure  15d.  Average  Power  Responses  f^r  Dolph-Chebyshev 
-30  dB  Weighting;  Planar  Array;  =  0.01,v  =  u/4 


Figure  15.  (Cont'd)  Average  Power  Responses  for  Dolph-Chebvshev 
-30  dB  Weighting;  Planar  Array 
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figure  15e.  Average  Power  Responses  for  Po 1 ph -i  heb\ she v 
-30  dB  Weighting;  Planar  Array;  v  =  u 


Figure  IS.  U.ontM)  Average  Power  Responses  for  Po l ph -t'hcby she 
-30  dB  Weighting;  Planar  Array 


Figure  17.  Average  Power  Responses  for  Hanning  Weighting;  Planar  Array 


Just  as  for  the  linear  array ,  a  simple  expression  for  the  saturation  level 
in  the  deep  side  lobe  region  may  be  obtained.  From  (25),  for  a  general  array 


(61) 


In  the  deep  side  lobe  region,  | A [  goes  to  zero  and  we  have 


?  = 


Y~  V 
*  K  vk 


k  x 


\r  r  v  M 

th  x  y£  9 


(62) 


giving 


e\?}  - 


k  X 


V  V* 

K I  VK  \ 


(65) 


For  uncorrelated  equal -variance  gain  perturbations,  (63)  simplifies  to 

£{f]  =  C  ?Whl  .  (641 
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Therefore,  for  w^  0,  the  saturation  level  for  the  mean  power  response, 
relative  to  the  peak  response  (for  0^  =  0),  is 


EfPl 

ic 


o' 


v<. 


<T2 


CfV  '  K.« 


(65) 


This  result  applies  to  an  arbitrary  array,  not  necessarily  planar. 
For  a  planar  array,  as  described  earlier, 


K, 


elf 


.W. 


(66) 


Y\ 

and  if  the  weights  are  multiplicative  (see  (19)),  then 

l^eff  ^  eff  eff 

- -  V*l 


(67) 


To  deduce  the  variance  of  P  in  the  deep  side  lobe  region,  we  use  (62) 
to  obtain 


V-V*V»V*  ' 


(68) 


For  Gaussian  which  are  uncorrelated  and  of  equal  variance, 


\\ 


rw<  ^  ^  +  +  ^8*3  ■ 


(69) 


When  this  result  is  employed  in  (68)  and  coupled  with  (64) ,  there 
follows 


Vorff]  =  +  |^v,‘|2 


(70) 


In  the  deep  side  lobes,  the  second  term  of  (70)  is  negligible,  and  we  have 


JRtev.ffl  *  <rr*  2lv,|‘  -  E{p], 
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Thus  CT  =  m,  giving  m  +  kf  =m(l+k),  and  (51)  and  table  2  hold  again.  These 
results  hold  for  arbitrary  arrays  with  any  weight  structure.  The  two  funda¬ 
mental  parameters  governing  behavior  in  the  deep  side  lobe  region  are  (T^ 
and  in  (65) . 


SUMMARY 

Results  for  the  mean  power  response  and  the  standard  deviation  of  the 
power  response  of  an  array  (linear  or  planar)  have  been  derived  for  arbitrary 
multiplicative  beamformer  weighting.  The  degree  of  widening  of  the  main  lobe 
response  caused  by  deep  side  lobe  designs  can  be  evaluated  quantitatively. 

A  simple  rule  of  thumb  for  the  saturation  level  in  the  deep  side  lobe  region 
has  been  presented;  it  depends  only  on  the  relative  variation  of  the  individual 
element  gains  and  on  the  effective  number  of  elements  in  the  array.  Specifica¬ 
tions  of  tolerances  of  element  gain  perturbations  can  be  made  quantitative 
once  the  array  size  and  type  of  array  pattern  have  been  selected. 
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Appendix  A 

PROBABILITY  DISTRIBUTION  OF  POWER  RESPONSE 

The  power  response  is  given  in  (25)  and  (7)  as 

T*  |^(l+0 'if  *  U+x  +  '.3f-(A>.«),4-(Ai*-j)')  {A.„ 

where 

*  +  "j  H  Y*V*‘  (A- 2) 

Now  x  and  y  are  correlated  zero-mean  real  Gaussian  random  variables.  In 
particular,  from  (A-2)  and  (24), 

e{(* + ^y]  *  f + 1 2  =  ^  S  \  vj  -  ^  > 

ej|«;3r}  =  \ v;  i  t>. 

(A-3) 

Therefore 

Eh'}  -  , 

Ety}  =  ±  &>-*)•  ■ 

EW-is.  5  ^V' 

1  (A-4) 

where  subscripts  r  and  i  denote  real  and  imaginary  parts,  respectively. 

2 

Then  the  probability  that  random  variable  P  is  less  than  T  is  given  by 
(A-l)  as 


A- 1 
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fv.V(?<  T‘)  =  ?r„l  ((A^xT  +  k*a)'  <  T‘) 


-  u\  [«'+»<  t)  =  (f  &  f  p(^-  (as 

*  <24^<TX 

where  the  probability  density  of  the  linearly  transformed  variables  o(-  Ar  + 
|3  =  +  y  is  joint  Gaussian: 


x, 


(A-6) 


Thus  (A-5)  corresponds  to  integration  of  a  bivariate  elliptical  normal  density 
over  a  circle  offset  from  its  peak  at  -  Ar,  A^.  After  rotating  coor¬ 

dinates  in  (A-5)  and  (A-6)  to  obtain  independent  random  variables,  we  can 
use  the  table  in  reference  1  to  evaluate  (A-5). 


A- 2 
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Appendix  B 

CHARACTERISTIC  FUNCTION  OF  POWER  RESPONSE 

From  (A-l)  through  (A-4) ,  we  have  the  power  response 
?  *  (Ar+x)’ 


and  the  probability  density 


The  characteristic  function  of  P  is  then 


-  (2-^ 9^1/^)  _  (B  3) 

where 

*  -  [2<('Y’)]  '  -  if 

f'  h’o  /)T’"? 

y*  6-1*5] 

jx  -  \  f  2  K 

p  =  i  2  A,  .  (B-4) 

Now  (reference  4,  3.323  2) 


B-l 


(B-5) 


^  ex^[*  4 / *  4  yy] 


1 

IT 

"1 

4^-r) 

,<>0  >fr*0,Kfr>V. 


All  the  conditions  in  (B-5)  are  satisfied  since  * ;  From  (A-3) f 

Isl  =  lEf(x+'bf]|  *  Eflx*i,r}  =  1>: 

therefore,  from  (A-4)  and  (B-6) , 

.  _  _ s<*  < , 

f  '  5*-  s;  ‘  S’+(D*-|S|') 

Employing  the  definitions  in  (B-4) ,  we  have 


i'X  4<r,'V  0  Y‘) 


and 


x  X  ('  X1  ^  A(  P  ‘  k  \* 

fr  **’  +«r"  *  W  +j  ar  - 


(B-6) 


(B-7) 


(B-8) 


(B-9) 


Combining  (B-3) ,  (B-5) ,  (B-8) ,  and  (B-9) ,  we  can  express  the  characteristic 
function  of  P  as 


f(f)= [i-if2(<+si) 


Cxy> 


I A l1-  i  f  2  (A*  n'+frfa1-  2  Ar  A,-  fi) 

I-  i  f  2 ^  <r; 


(B-10) 
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Now,  by  eliminating  variables  in  favor  of  the  original  quantities  via 
(A-4)  ,  there  follows 

2K+0-  2t> 


4  <  ri‘  i'-f) 


(B-ll) 


and 


-r'A 

f [f)  =  jt-  if 2 D  *(j ^ (tf-lsfljj  evj, 


if 


1-  if2P  +(i'f),(t)2-]Sr) 


CB- 12) 


A  single  numerical  Fourier  transform  in  (B-12)  would  yield  the  probability 
distribution  of  ?  (reference  2)  for  specified  values  of  A,  D,  and  S.  We  do 
not  pursue  that  tack  but,  instead,  derive  the  cumulants  of  P.  In  (B-12),  we 
let 


X=if}a  =  2P,  \>  -  DMs)\  c=)M\  I Af  I>-  ?eKs*], 

where  a>0,  b>0,  c>0,  d  >  0,  to  obtain 


r  ^  >T,/4  1 

"  C-x<i 

h-  XQ4-X  t>j  txf 

X  *1 

l-x«+xl>J 

We  then  expand  J0v\f  in  a  power  series  in  x, 

x"' 

*  *  ■  ’ 
obtaining  the  cumulants 


CB— 13) 


(B- 14) 


C  B- IS) 


(B-16) 


B-3 


where 


M 

k*o 


(rD  4h-k\^-2kLK  ^  , 

i«_k  \  w  J  ^  b  ,  ^ 


.  k 

l  , 


CB— 17) 


In  particular, 

c,  «  2  a  +  c 

r  =  J-cf-  \>  +  2ca-2d 

'•j  i  ^ 

c  *  cx  —  +  (o da  . 

3  (B-18) 


In  terms  of  the  original  quantities,  using  (B-13) , 

C,  *  IA\*+  T>  -  Efp] 

c,*  t?  4  |sf42D|M’  +  ^A,S*}  *  Vorf ^ 

q.2p’+CI>lSl'+£p2lAf4  C|S|’|A|‘+  ODRA’S*} 

(B- 19) 


The  first  relationship  in  (B-19)  for  the  mean  C2  of  P  actually  holds 
for  any  distribution  of  {r^},  not  just  Gaussian;  this  is  most  easily  seen  by 
direct  use  of  (25),  (24),  (7),  and  (A-3) . 
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MOMENTS  OF  POWER  RESPONSE  FOR  NON-GAUSSIAN  PERTURBATIONS 


From  (25)  and  (7),  the  power  response  is 

?  =  |A  +  f  r*v*|\  (C-D 

The  mean  of  P  was  derived  in  (32)  and  is  given  by 

=UI‘  ■*->»,  I1  .  (c-2) 

The  mean  square  value  of  P  is  available  from  (C-l),  (31),  and  the  independence 
of  {r^},  as  follows: 

El?']  =  A  * f  *  JS v/)(A -vSc v*)| 

*  lAl  +  Af  +  A1  **  +  f\*  K 

+2 +2A  ^ +  ^{'kv« r^’r"]^Kyt • 

(C-3) 

The  fourth-order  average  in  (C-3)  has  the  value  m4  if  all  four  subscripts 
are  identical,  and  the  value  if  two  of  the  subscripts  are  identical  and  the 
other  two  are  identical,  but  not  with  the  former  pair.  Then  the  last  term 
in  (C-3)  is  expressible  as 

wi,  f v*'4  +  <  [?  I  w  m  *  %■<  v/] 

U  k**  k*f  J 


(C-4) 


Equation  (C-3)  then  becomes 

dH  -  lAlV-tlAlVf  M’  +  2«,1?e?Vf '?’) 

+  •4wi3Re[A2|\(,l’v,*]  v  «4  Xlv.1* 

Combining  this  result  with  (C-2) ,  there  follows 

Vovff]  »  2W,[|Arflv.|*  +  Ke|A*^}] 

+  4»,K*lA  Zhrfv,*] 

+  [(S-KI*)'  + 1 f]  +  K-3<)XKf. 
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Appendix  D 


ELEMENT  GAIN  PERTURBATIONS  EXPRESSED  IN  DECIBELS 


The  voltage  gain  of  the  k-th  element  was  given  in  (23),  (7),  and  (8) 


as 


(l  +  *iK  *  e» 


(D-l) 


Then  (for  w^  >  0)  the  probability 

TU(e„<V)  =  ?rA(n(l+0<V)-f,»V(r-<'|-')'-T«(^'l),(D'2) 


where  is  the  cumulative  distribution  function  of  random  variable  r^. 


Choose  voltage  levels  V^,  such  that 

-  X  ('fti  '■■>), 

PrA(=s<VL)  -  rL  (*ai  «•»•). 


That  is. 


T  ( =  T 
'k  V  y  U)L  1 


or  voltage  levels 


V 


CD-3) 


CD— 4) 


(D-5) 


These  are  the  voltage  levels  within  which  element  voltage  gain  e^  can  be 
found  with  tail  probabilities  T^,  1  -  T^. 

We  now  define  decibel  levels  (corresponding  to  the  specified  probabil¬ 
ity  values  Ty  ,  relative  to  the  design  value  w^,  as 


D-l 


IR 


The  range  (JBj  ,  dByj  is  that  range  within  wh  i  v  h  the  element  gain  per 
tions  can  be  found  with  tail  probabilities  I|,  1  !y.  i  Ho  get  tin 

result  if  we  consider  the  power  gain  ej  as  the  fundament.il  r.uuion 
provided  that  — 1,  w^  *  0.)  If  ^arlr^f  -  o,  then 

S'  *  ?  <\ 
l  0,  v.  o' 

X  (uW  0  ft*  0<  JP  I, 

and  dB^j  ^  =  0,  as  expected. 

For  Gaussian  element  perturbations, 


where  £  is  the  Gaussian  cumulative  distribution  function 


$  (f)  -  {  (?-*-)  e  Xf  (-  uVi) 

— 

Then  (D-6)  yields 

dB^-  20  loj(l+«J 

where  $  is  the  inverse  $  function. 

For  uniform  element  perturbat ions* 


T,W*  *fST  +  i  ’  Ixl</Ts; 

ifW)-  2t/3V,  (3- ,  0<y<  I 


D-2 


t  !J  If  M 

‘-arm* 

V  a  r  1  af’  1 


n-  s 


I  I)  9 


I  b-  I  t» 


ID-  1  1 
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I  hen  ;  L) -  M  v  i  e  1  d s 


CD-12) 


filiations  !  M-  1 0 1  and  (D-IJ)  are  plotted  in  figures  3  and  4,  respectively. 
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